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ABSTRACT 

The tight relation of star formation with molecular gas indicated by observations and assumed in 
recent models implies that the efficiency with which galaxies convert their gas into stars depends on 
gas metallicity. This is because abundance of molecular hydrogen is sensitive to abundance of dust, 
which catalyses formation of H2 and helps to shield it from dissociating radiation. In this study 
we point out that in the absence of significant pre-enrichment by Population III stars forming out 
of zero metallicity gas, such H2-based star formation is expected to leave an imprint in the form of 
bi-modality in the metallicity distribution among dwarf galaxies and in the metallicity distribution 
of stars within individual galaxies. The bi-modality arises because when gas metallicity (and dust 
abundance) is low, formation of molecular gas is inefficient, the gas consumption time scale is long, and 
star formation and metal enrichment proceed slowly. When metallicity reaches a critical threshold 
value star formation and enrichment accelerate, which leads to rapid increase in both stellar mass 
and metallicity of galaxies. We demonstrate this process both using a simple analytical model and 
full cosmological simulations. In contrast, observed metallicity distribution of dwarf galaxies or stars 
within them are not bi-modal. We argue that this discrepancy points to substantial early stochastic 
pre-enrichment by population III stars to levels Z ~ \Q~ 2 Zq in dense, star forming regions of early 
galaxies. 

Subject headings: cosmology: theory - galaxies: evolution - galaxies: formation - stars: formation - 
methods: numerical 



1. INTRODUCTION 

The observed correlation between stellar mass and 
metallicity, established for a wide range of galaxy masses, 
morphologies, and metallicities, is an important diag- 
nostic in the investigation of physical processes that 
shape the baryonic component of galaxies. iDekel fc Wool 
(2003) find a correlation of the form Z oc M" z with 
olz ~ 0.4 for the Local Group d warf spheroidal galax- 
ies. Anal ysis of the SPSS data , by lTremonti et all (|2004f ) 
(see also iGallazzi et al.l 120051) indicates that the slope 
of the scaling increases with decreasing M*: it is quite 
flat for M» > 3 x 10 10 M(T), while it st e epens for lower 
masses and exceed s the IDe kel & Woo (2003) value by 
M ~ few x I0 8 M Q . iLee et al.l (|2006l L on the other hand, 
find that at even lower stellar masses the scaling becomes 
less steep, with a slope of az = 0.29 ± 0.03 for the dwarf 
galaxies in the ir sample ( 10 6 M fi < M* < few x 10 8 M ). 
Most recently. iKir bv et al.l (|201lD find that at even lower 
stellar masses (M* < 10 5 M Q ) the relation becomes flat, 
a z ~ 0. 

The stellar mass-metallicity relation of low-metallicity 
systems is sensitive to details of star formation, en- 
richment, and metal transport during the earliest 
stages of galaxy evolution. For example, it was argued 
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that Af* — Z is sen sitive to blow-out of metals via 
galac t ic winds (e.g.. ITassis et all 120031 : IDekel fc Wool 
20031 IKobavashi et al.l I2007D or to increasing ineffi- 
ciency of star form a tion with decr e asing galaxy mass 
dBrooks et al.l l2007t Ide Rossi et all 120071: ITassis et all 
120081: iGnedin fc Kravtsovl l2011aft The latter can be 
due to star form ation regulation by feedback (e.g., 
iBrooks et al.l l2007) or by processes controlli ng formation 
of de nse, star forming regions within ISM (Tassi s~et all 
2008]). One such key process is regulation of molecular 
;as formation by dust abundance. Both observations 
Wong fc Bllti l'2002at iLerov et al.l [20081 iBigiel et all 



120081: iGenzel et all l2010UBigiel et all 120111) and the 
oretical considerations (e.g., iKrumholz "etal! f20Tlh 
suggest that star formation is tightly correlated with the 
mass of gas in molecular phase. The latter is sensitive 
to dust abundance, which catalyses the formation of 
H2 an d helps to shield i t from d issociating radiation 
(e.g., iSpaans fc Norman] 119971 : iHirashita fc Ferraral 
2002: IKrumholz et al.l 120091: IGnedin et al l l200a 



Gnedin fc Kravtsovll2011bl lAvkutalp fc Spaansll201lir 

One of the implication of this tight relation is a "catch- 
22" situation that exists during early stages of galaxy 
evolution and chemical enrichment: stars are needed to 
produce metals and dust, while dust is needed to facil- 
itate star formation. Star formation can thus be ex- 
pected to proceed first in the central regions of galax- 
ies, where gas surface density is high and it is easier 
for H2 to form even with small amount of dust. As 
central regions are enriched, the metals can be mixed 
with gas at larger radii either by internal mixing pro- 
cesses or via mergers of gaseous disks. As lower sur- 
face density regions are enriched, formation of molecu- 
lar gas and star formation will become more efficient at 
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larger radii within galaxies. Given that this sequence 
requires a certain amount of time, which can be pro- 
longed by loss of metals via feedback-driven winds, star 
formation in galaxies can be expected to be delayed or 
even completely suppressed. The delay or suppression 
depend on the surface density of gas within galaxies and 
can thus be expected to be the largest in dwarf galax- 
ies, whic h have the smallest gas surface densities. Sim- 
ulations (|Gnedin fc Kravtsovll20"llat iKuhlen et al.ll20H"l) 
and analytical models iKrumholz fc Dekell (|2011l ) show 
that strong suppression of star formation via such pro- 
cesses occurs in halos of total mass Mh < 10 10 M . 

In this paper we examine how such processes affect the 
metallicities of galaxies and the metallicity distribution 
of stars within them in more detail by using both simple 
analytical model (§ [2]) and full cosmological simulations 
(§[3]) based on the model for H2 formation and H2- based 
star formation introduced in iGnedin et al.1 (|2009f ) and 
iGnedin fc Kravtsovl (|2011bf ). We show that when star 
formation is tied to molecular gas, bi-modality in metal- 
licity arises generically due to the existence of a charac- 
teristic gas surface density below which star formation 
efficiency is suppressed. However, as we discuss in § [4J 
observations do not show evidence for such bi-modality 
in high redshift or local low-metallicity systems. We in- 
terpret this apparent discrepancy between theoretical ex- 
pectation and observations as indication that gas within 
regions where first Population II stars were formed was 
pre-enriched by Population III stars to metallicities as 
high as ~ Zq. Such interpretation is consistent with 
models of dwarf Milky Way satellites, the metallicity dis- 
tribution of the Milky Way halo stars, and observations 
of metallicity distribution of the damped Lya absorp- 
tion systems, which also indicate th at pre-enrichment 
at least to Z ~ 10~ 3 Zq is required (iHelmi et al J 12006b 
Rollinde et all 120091 : iKirbv et all 120111 : iKarlsson et alJ 



201 If) . Furthermore, abundance patterns in very low 



metallicity systems are in good agree ment with those 
expected from the Popu l ation III stars (iBromm fc Loebl 
[2001 ISimon et alJl2010b iGooke et all 1201 lfi . Significant 
pre-enrichment is also suggested by simulations of feed- 
back by the Population III stars ()Wise et alJl201lh . 

2. METALLICITY-MEDIATED STAR FORMATION AND 
ENRICHMENT HISTORY OF GALAXIES 

We begin by illustrating how bi-modal metallicity dis- 
tribution can be expected to arise generically both for 
the distribution of metallicities within galaxies and the 
distribution of average metallicities of galaxies. We do 
so using a simple analytical model which, nevertheless, 
captures the relevant physics in a qualitative way. 

Studying the env ironmental effects o f H^based star 
formation models, IGnedin fc Kravtsovl (|2011aD found 
that the surface density of the star formation rate as a 
function of the total neutral gas surface density, Eh i+h 2 > 
can be well approximated as 
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where Ehi+h 2 is the neutral gas consumption rate per 
unit area and E cr ; t is a characteristic "threshold" surface 
density of neutral gas below which the Esfr — Ehi+h 2 
relation steepens as gas becomes predominantly atomic 



at lower surface densities. Defining E t = Ehi+h 2 + E*, 
where E* is stellar surface density, the (neutral) gas frac- 
tion at each location is f g = Eni+H 2 /Et, and the evo- 
lution of f g at constant E t is described by the following 
differential equation: 
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In general, the threshold surface density E cr ;t is a com- 
plex function of the dust-to-gas ratio (which we assume 
to be linearly proportional to the gas metallicity Z in 
this paper), and the interstellar radiation field. However, 
for low Z, E crit is approximately inversely proportional 
to Z ', with the coefficient of proportionality being some- 
what dependent on the value of the interstellar r adiation 
field [see Figure 6 of lGnedin fc Kravtsovl (|2011al) ]. Radi- 
ation fields in high-red shift galaxies are g enerally higher 
than in local galaxies ([Chen et all 12009 !). so in the fol- 
lowing we assume E cr jt ~ lOM0pc~ 2 /Z, which is ap- 
propriate for radiation fields 10 to 100 times higher than 
in the Milky Way. Since E cr i t depends only on approxi- 
m ately the square root of the radiation field [see eq. 14 
of lGnedin fc Kravtsovl (]2011al )]. the range above corre- 
sponds to a factor of 3 uncertainty in E cr j t . 

The gas fraction and the gas metallicity are related 
through 

Z = Z + t/ e ffln(l//g), 

where Zq is the "primordial" metallicity (i.e. the metal- 
licity from star formation that is uncorrelated with the 
dust abundance) , whic h may reach values exceeding 10 -3 
solar (|Wise et al.ll2011l ) , and y e ff is the so-called "effective 
yield" , which, in a closed-box approximation, is equal to 
the true metal yield. Assuming that the effective yield is 
a weak function of time (or, equivalently, gas metallicity), 
we can derive a single equation for the gas metallicity as 
a function of time, 
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where q = E t /(10 M© pc~ 2 ). 

The solution to Equation ([3]) can be obtained analyti- 
cally for the limit Z <C y e e, if we introduce a new time 
variable dx = (y c s / T SF)dt, which "hides" all complex de- 
pendence of tsf and y c s on other physical quantities. In 
general, however, tsf will depend on the gas metallicity, 
so finding a solution of Z as a function of x only produces 
an implicit dependence of Z on the actual physical time 
t. For simple esti mates, however, we can assume tha t 
t sf w 1.5-2Gyr (jBigiel et alll201lHGenzel et al.|[2010h . 

With these simplifying assumptions, Equation <j3j) re- 
duces to 

1 V 2 



dx 



whose implicit solution is 
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Solution ([S]) exhibits different behavior for different 
values of Zq and q. If qZQ -c 1, the metallicity re- 
mains just above the "primordial" value Zq for a long 
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Fig. 1. — Top panel: metallicity (in solar units) as a function of 
the total surface density at At = rgp. Bottom panel: time interval 
required to reach a particular metallicity as a function of the total 
neutral gas surface density. The val ue of y e ff = Q.2Z p) is assumed 
in the presented model calculations (Dalcanton 2007). 

time, x ^> 1, keeping galaxies in the low-metallicity, 
low-stellar-mass part of the parameter space. However, 
if qZq 3> 1, Z is growing approximately linearly with 
time, Z w x + Zq, and the gas metallicity reaches the 
value Z = y c ff (recall that the solution is only valid 
for Z < y c ff) in a time interval At = tsf- Finally, for 
Z > y c g, the exponential term in Equation ([3]) becomes 
non-trivial, and the growth rate of Z slows down consid- 
erably and eventually becomes only logarithmic. Figure 
Q] illustrates the solution J5} as a function of the total 
surface density for Z at constant t and t at constant 20. 
The transition between the two regimes (low E t , low Z 
and high S t , high Z) is rather rapid when the value of 
Zo is low, so the two regimes are well separated. The 
transition becomes shallower as Zq increases. 

It is this transition that is the source of bi-modality in 
the distribution of galaxy metallicities when Zq is low. 
Given that surface density of gaseous disks is expecte d 
to decrease with decreasing halo mass (|Mo et alj |1998). 
halos below a certain mass will host disks with surface 
density smaller than the above transition value. Star 
formation and metal enrichment would proceed in such 
disks at a slow rate and they would stagnate in the low 
M*, low Z regime. Conversely, gas disks with high sur- 
face densities, hosted by sufficiently massive halos, will 
form substantial mass of stars and will reach high metal- 
licity within period of < rgp « 1.5 — 2 Gyrs, at least 
within their central regions. A qualitatively similar ef- 
fect can be expected for the metallicity distribution of 
stars within individual galaxies, for systems that spend 
significant fraction of their evolution in the stagnant, low- 
metallicity phase. While a galaxy is in this phase, metal- 
licity evolves slowly and stellar mass that is building up 
will predomninantly consist of low Z stars. If at some 
point in its evolution either the surface density of its gas 

6 Changing £ cr it in our solution would translate the curves hori- 
zontally by the corresponding factor, but it would not change their 
shape. 
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Fig. 2. — Mass-metallicity relation for simulated galaxies at dif- 
ferent redshifts (symbols/colors as in legend). 

disk or the metallicity increase to the level required by 
the condition above, the galaxy will undergo dramatic in- 
crease in star formation rate and will rapidly self-enrich 
itself producing a population of high Z stars. Given that 
the transition is rapid, there will be a relative paucity 
of stars at intermediate metallicities and the metallicity 
distribution will be bi-modal. Note that pronounced bi- 
modality can be expected only for galaxies that produced 
significant fraction of stars in both regimes. Clearly, if a 
galaxy has always evolved in the stagnant regime, it will 
only have low metallicity stars without a high Z peak, 
while the low Z peak may be insignificant in a galaxy 
that formed most of its stars after the transition. 

The simple model and arguments presented above 
show that bi-modality in metallicity distribution is gener- 
ically expected in galaxies of a certain mass range. We 
now turn to simulations that include the process of disk 
formation and H2-based star formation self-consistenly 
and demonstrate that the expected bi-modality does in- 
deed arise naturally for simulated galaxies. 

3. METALLICITY BI-MODALITY IN SIMULATIONS WITH 
H 2 BASED STAR FORMATION 

To examine detailed predictions for the metallicity dis- 
tribution, we use a set of high-resolution sim ulations 
of WMAP$] cosmology (jKomatsu et all 12009ft carried 
out with the Adaptive Refinement Tree (ART) cod e 
(|Kravtsovill999l: iKravtsov et "al]|2003 Hudd et al.ll2008ft . 
The detail ed description of the simulation is presented 
elsewhere dZemp et al.H20li[) . The simulation volume of 
25.6/i _1 Mpc includes several regions of interest which 
are sampled in the initial conditions with the effective 
resolution of 2048 3 , resulting in the dark matter particle 
mass of 2 x 10 5 M Q . Each region of interest is allowed to 
refine dynamically up to the comoving spatial resolution 
of 200h~ 1 pc, which is equivalent to the spatial resolution 
in physical units of (57, 70, 90, 95) pc at z = (4, 3, 2.5, 2). 
Regions of interest are chosen from a lower resolution dis- 
sipationless simulation to include 5i? v ir spheres around 

7 fi M = 0.28, U B = 0.046, h = 0.7, n s = 0.96, ag = 0.82. 
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Fig. 3. — Metallicity distributions of stars in three simulated galaxies with different stellar masses. The distributions in the larger 
and medium-mass galaxies (left and middle panels) are clearly bimodal.The extremely-low stellar mass galaxy does not have any stars of 
metallicity higher than 10 — 3 solar. 



six galaxies with z — masses between 10 12 and 10 13 
solar masses. 

The physi cal processes modelled in the simu lation are 
described in IGnedin k Kravtsovl f2010t 12011 rl . In par- 
ticular, we include gas cooling, including cooling by met- 
als, molecular hydrogen, and dust, a phenomenologi- 
cal model for molecular hydrogen formation, full time- 
dependent and spatially variable 3D radiative transfer 
of ionizing and Lyman- Werner band radiation using the 
Optically Thin Variable Eddingt o n Ten sor (OTVET) ap- 
proximation of iGnedin k Abel ) (120011) . star formation 



in the molecular gas using the iKrumholz k Tanl (f2007T l 
recipe, and (in our fiducial run) the supernova feedback 
modelled as thermal energy injection over 20Myr after 
the formation of a stellar population. 

Our simulation does not resolve the formation of the 
Population III stars and associated enrichment. How- 
ever, this initial enrichment and the dust that results 
from it are required for a more widely spread formation 
of molecular gas and the Population II star formation. 
Therefore, we model the enrichment by the unresolved 
Population III stars by assuming that a gas dense enough 
to form molecules has a minimum amount of cosmic dust 
corresponding to dust-to-gas ratio 10 -3 times smaller 
than the Milky Way value (i.e . , Dmw = 10~ 3 in the 
notation of IGnedin k Kravtsovl 12011 bp ). Given that we 
assume that the dust-to-gas ratio scales linearly with gas 
metallicity, this assumption is equivalent to setting the 
initial metallicity to Z$ = 1O _3 Z0 in the notation of the 
previous section. This "dust floor" is only accounted for 
in computing the molecular hydrogen abundance. It is 
not added to the actual metal abundance followed in the 
simulation and is not taken into account in calculation 
of cooling rates, for example. However, in the discus- 
sion below we will examine results by introducing explicit 
metallicity floor during post-processing of simulations. 

Star formation in the simulations is b ased on abun- 
dance of m olecular gas, as des c ribed i n IGnedin et al.l 
(|2009( l and IGnedin k Kravtsovl (|2011bf ). Our key re- 
sult crucially depend on the assumption that star for- 
mation is tightly correlated with H2 abundance. The 
conclusions can therefore be subject to uncertainties as- 
sociated with uncertainties of H2 modelling. Note that 
one of the key ingredients of the H2 formation model is 



H2 formation rate on dust, for which we adopt the em- 
pirically constrained rate of iWolfire et al.l (|2008| ) and do 
no t assume any temperature d ependence (see Appendix 
in IGnedin k Kravtsovl 1201 lbl for details) . It was ar- 
gued, however, on theoretical grounds that such temper- 
ature dependence may exist ([Burke k Hollenbachl 119831 : 
ICazaux k Spaans 2009). Both the degree of the temper- 
ature dependence and its effects on our conclusion are 
uncertain and such uncertainties of molecular hydrogen 
formation in environments drastically different from the 
Milky Way should be viewed as one of the caveats to our 
results and conclusions. 

The mass-metallicity relation for the model galaxies 
from our fiducial simulation is shown in Figure [21 In 
concord with the expectations presented in fJ21 the dis- 
tribution of simulated galaxies is clearly bimodal: galax- 
ies with stellar masses > few x 10 8 M Q (residing in halos 
with total masses larger than about ~2x 10 10 M Q ) and 
metallicities higher than w 0.1 Zq comprise a distinctly 
different population, with a much weaker scaling between 
stellar mass and metallicity compared to galaxies with 
stellar masses < 1O 7 M and metallicities < 10~ 2 Zq. 
Galaxies with intermediate values of stellar mass and 
metallicity are scarce, but they appear to be following 
the well-defined extension of the M* — Z scaling of lower 
mass galaxies. Interestingly, the results do not change 
significantly with decreasing redshift: the low-M*, low- 
Z population does not become significantly depleted with 
increasing cosmic time. 

Note that this bi-modality in the distribution of galaxy 
metallicities is a unique signature of H2-based star for- 
mation, set by the sensitivity of H2 abundance to the 
abundance of dust. Such bi-modality is not present in 
simulations which use common star formation recipes 
based on the total gas density rather than o n the density 
of on ly molecular gas (see, e.g., Fig. 6 in iTassis et al.1 
2008). Dust is important both as a catalyst for H2 for- 
mation and as a shield from dissoc i ating UV radiation. 
As shown by lAvkutalp k Spaansl (|201lD , UV heating 
may also suppress formation of dense regions in the first 
place. Our results indicate that such suppression in high- 
rcdshift small halos is not complete. Overall, however, 
such heating contributes to increasing inefficiency of gas 
conversion into stars in such systems and thus is one of 



the reasons for the bi-modality that we observe. 

Figure[3]demonstrates that bi-modality can also be ap- 
parent in the metallicity distribution of stars within in- 
dividual galaxies. It shows the metallicity distribution of 
stellar particles in three galaxies of different stellar mass 
at z — 3. The galaxy in the left panel has a stellar mass 
of 3 x 10 10 M©, and its distribution of stellar metallicities 
is bi-modal, although the fraction of stars in the low- 
metallicity peak is small. The distribution of metallici- 
ties in galaxy of intermediate stellar mass (1.7 x 10 8 M Q ), 
shown in the middle panel, exhibits distinct bimodality 
with the low-Z peak containing a significant fraction of 
stars. Finally, the smallest galaxy with stellar mass of 
1.5 x 10 6 M©, shown in the right panel, has only a narrow 
\aw-Z peak. This galaxy never evolves to undergo a tran- 
sition to the fast star formation regime and thus never de- 
velopes the high-metallicity peak present in larger mass 
galaxies (see discussion in the previous section). 

4. COMPARISON WITH OBSERVATIONS 

Color points in Figure [4] show the M» — Z correlation 
for the Local Gro up dwarf galaxies (jKalirai et al.l l2010t 
iKirbv et al.1 1 2 1 lh and other nearby dwarfs fLee et all 
2006) covering a mass range of M* » 10 3 — 10 8 M Q 
alon g with the mean rela tion for more massive galax- 
ies dTremonti et al.l 120041 ) . The metallicities are mea- 



sured for the Local Group dw arfs using stellar iron abun- 
dance, while iLee etaLl (|2006f ) and lTremonti et al.l (|2004f ) 
measured gas oxygen abundances. We rescale th e latter 
using the empirical prescription of iSimon et al.l (2006) 
[eq.(12)]. After such rescaling both the Local Group 
dwarf spheroida ls and dwarf irregular galaxies in the 
ILee et afl (|2006l ) sample exhibit a fairly tight correla- 
tion that joins smoothly into the relatio n measured for 
the la rge sample of the SDSS galaxies bv lTremonti et all 
(|2004f ). Note that Z — relation flattens substantially 
at the lowest stellar masses: at M* < 10 5 M Q metal- 
licities no longer exhibit a tight correlation with stel- 
lar mass and are instead scattered in the range Z ~ 

10~ 2 - 5 - io- 1 - 5 . 

Even a cursory comparison of the predicted relation 
shown in Figure [2] with the correlation of real galax- 
ies in Figure 0] clearly shows that the two are inconsis- 
tent. The main difference is that simulated dwarf galax- 
ies have drastically lower metallicities compared to the 
observed ones at a given stellar mass. This discrepancy 
can have different explanations and we discuss some of 
them below. Here we note that one possible explanation 
is that gas was universally pre-enriched to a metal floor of 
Z m in ~ 1 X 10~ 2 . Such pre-enrichment could indeed ex- 
plain the apparent lack of correlation between metallicity 
and stellar mass for the smallest dwarf galaxies. Given 
that different systems and regions in the early universe 
evolve somewhat differently, we can generally expect such 
pre-enrichment to be stochastic. Filled black circles in 
Figure H] show simulated galaxies in which a stochastic 
metallicity floor was added, which we model as a log- 
normal with median metallicity of O.OIZq and scatter of 
&iog w z = 0.2 dcx. Although the agreement of observa- 
tions with such floor is not perfect, it is much improved 
compared to when the original simulation metallicities 
are used. Substantial pre-enrichment of gas of dwarf 
galaxy disks by preceding Pop III stars can thus help 
explain the discrepancy between predictions in Figure [2] 
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Fig. 4. — Stellar mass-metallicity relation for nearby dwarf galax- 
ies. Blue crosses: iron metal licity observation s in local dwar f galax- 
ies from Kirby ct al. (2011); magenta stars: Kalirai ct al. (2010); 
orange stars: dwarf galaxies from sample of ILee et al.l 120061). Av- 
erage relation measured for the SDSS galaxies by Tremonti ct al. 
(2004) is shown by colored dashed line. Simulation results from 
2 = 3, with added stochastic metallicity floor described by log- 
normal pdf with the median of 1 X 10 — 2 solar and scatter of 0.2 dex 
are shown for comparison with the black solid dots. 

and observations in Figure |4j 

The distribution of metallicities of dwarf galaxies in 
Figure [4] also does not exhibit any clear bi-modality. 
Observations do not seem to show a paucity of galax- 
ies at stellar masses ~ 10 6 — 10 8 M Q and metal- 
licities Z ~ 0.01 — 0.1. The interpretation of this 
discrepancy is less certain, given that the sample of 
dwarf galaxies shown in the figure is not volume lim- 
ited. Nevertheless, the lack of bi-modality in the distri- 
bution of average metallicities of dwarf galaxies is con- 
sistent with the lack of corresponding bi-modality in the 
metallicity distribution of stars within galaxies. These 
distributions have been measured for both the Milky 
Way and for man y of the dwarf spheroidal galaxies in 
the Lo cal Group |Beers fc Christliebl 120051: iHelmi et all 
20061: iStarkenburg et al.l 120101 iBattaglia et al.1 120111: 
Kirbv et al.l l2011t iCarrera et al.l |2011|) with "the statis- 
tical power that would be sufficient to uncover the level 
of bi-modality seen in our simulations. However, none 
has been detected so far. 

We elaborate on the possible causes of such a stark 
disagreement between the generic prediction of the 
molecular-gas-based star formation model and observa- 
tions in the next section. 

5. DISCUSSION 

As we show in §[3 and [31 the assumption that star for- 
mation is tied to the molecular phase of the interstellar 
medium naturally leads to bi-modality in the distribu- 
tion of stellar masses and metallicities of dwarf galaxies. 
Such bi-modality is not observed in the available data for 
nearby dwarf galaxies and in metallicity distribution of 
stars within such galaxies. 

Given the absence of bi-modality, which we argue is 
a generic feature of H2-based star formation, one can 
ask whether the assumption about tight relation be- 
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tween star formation and molecular gas is valid. Relax- 
ing such assumption could perhaps modify the predic- 
tion about pronounced bi-modality. Although in prin- 
ciple other models that do not consider molecular gas 
as a primary indicat or of star forming gas are possible 
(jOstriker et al.l [20101 e.g.), we note that there is cur- 
rently good motivation for H2 based star formation be- 
cause star formation rate on ~ 100pc-l kpc scales corre- 
lates strongly with Ho mass in the c orresponding region 
(I Wong fc Blitz! I2002M I Bigiel et al.l [2001 IGenzel et al.1 
l201CtlBigiel et al1 l2011). Such star formation model also 
naturally explains the observed tre nd of the Kennicutt- 
Schmidt relation with metallicitv (|Gnedin fc Kravtsovl 
120101 IRafelski et al.H201lt IBolatto et al.ll2011fl . 

An alternative explanation for the observed lack of the 
bi-modality is primordial pre-enrichment of the ISM and 
part of the IGM by the first Population III stars and, 
perhaps, by rapid subsequent formation of Population II 
stars. Such pre-enrichment is often dubbed "the metal- 
licity floor" , although the observed distribution of metal- 
licities in the lowest mass dwarf spheroidals in the Lo- 
cal Group is better described by a lognormal distribition 
with the median at Z w 10 _2 Zq and an rms scatter of 
0.2 dex. Indeed, theoretically it is more natural to expect 
some probability distribution density of initial metallici- 
ties, given that different regions collapsing into a galaxy 
may have had different evolution ary histories and e nrich- 
ment levels (see, e.g., Fig. 2 in iWise et al.l 120111 ). We 
can thus expect metallicities to be somewhat different in 
different galaxies and in different regions within a given 
galaxy. 

We show that assuming a minimum metallicity 
stochastically varying according to the log-normal pdf 
with parameters given above significantly improves the 
agreement of simulated galaxies with observed low-mass 
dwarf galaxies. Although pre-enrichment to Z ~ 10~ 2 Zq 
may seem extreme, it is actually in line with current 
theoretical model s. For example, a recent study by 
IWise et "all (|2011[) shows that Population III stars can 
enrich gas of the entire host halo to metallicity of Z ~ 
10 _3 ^q. It is important, however, to distinguish mass- 
weighted average metallicity of all the pre-enriched gas 
and metallicity of the central regions that will be the sites 
of Population II star formation. The latter can be ex- 
pected to be enriched to metallicities c onsiderably higher 
than average. Indeed, simulations of IWise et al.l (|2011l ) 
show that although gas in dwarf halos at z > 10 is en- 
riched by Population III stars to average mass-weighted 
metallicities of ~ 10~ 3 — 10~ 25 Z Q , the first Population 
II stars may form in the central regions of halos with 
metallicities of Z* w 10~ 2 Z & or even higher (see their 
Figure 3). 

The pre-enrichment hypothesis can explain the sharp 
decline in the distribution of Milky Way halo stars be- 
low [Fc/H] ~ —3, the flattening of the mass-metallicity 
relation to a plateau at low metallicities, as well as the 
absence of bi-modal behavior in metallicities and other 
observed properties of the baryonic components of dwarf 
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galaxies. Indeed, if the gas in dense, star forming regions 
of early galaxies becomes pre-enriched by the first popu- 
lation of stars to a level of ~ 10~ 2 Z & the bi-modality of 
the galaxy metallicity distribution will become consider- 
ably less pronounced (see the blue line in the top panel 
of Fig. Q]). 

A related observational fact is the lack of extremely low 
metallicity (Z < 1O _3 Z0) Dampe d Ly-q (DLA) systems 
in current observat ional samples ()Prochaska et al.1 [20031 : 
iKaplan et all [2010. Rafelski, M. et al. 2011, in prepa- 
ration). Observations are sufficiently sensitive to detect 
systems with metallicities as low as « 10~ 4 Z Q , yet no 
systems with metallicity < 3 x 10~ 2 Zq have been yet 
found even in the highest redshift systems (4 < z < 5). 
The lack of low-metallicity DLA systems also points to- 
wards significant early pre-enrichment of galactic disks. 
It remains to be seen whether simulations modeling pre- 
enrichment by Population III and II stars will be able to 
reproduce such fairly universal level of pre-enrichment in 
different stellar systems and gaseous disks (e.g., no place 
in the universe is clearly demonstrated to have metallic- 
ity much below 10~ 3 solar). 

The hypothesis that most of the ultra-faint Local 
Group dwarfs (M* < 1O 6 M ) are pre-enriched by the 
first generation of stars implies that essentially all heavy 
elements observed in systems with M* < 10 5 M Q have 
been produced by the Population III stars. This argu- 
ment adds support to independent arguments for signif- 
icant Population III pre-enrichment based on chemica l 
evolution modelling of dw arf galaxies (iKirbv et al.ll2011l) 
and the Milky Way halo (|Rollinde et al.l 12009( 1. and ex- 
aminati on of abundance p atterns in stars belonging to 
dwarfs (|Simqn et al.ll201(J ) and in damped Lyq systems 
(|Cooke et al.ll201lL see also lKarlsson et aill2011l for a re- 
cent review). Measurements of the abundance patterns 
in spectra of stars within these galaxies can thus directly 
probe the nucleosynthetic yields of the very first super- 
novae. This opens an exciting avenue of using these 
smallest galaxies as unique window into the properties 
of first stars and conditions for star formation during the 
earliest stages of galaxy evolution. 
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